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EXECUTIVE  SUMMARY 


The  advent  of  imaging  seekers  in  the  3-5  pm  and  8-12  pm  bands  poses  a  new  challenge  in  the 
area  of  infrared  (IR)  countermeasures.  Area-extensive  infrared  countermeasures  will  be  necessaiy  to  defeat 
these  new  seekers.  Based  on  work  by  the  Army,  one  effective  component  of  these  area-extensive  IR 
countermeasures  may  be  obscurants.  These  materials  can  degrade  the  scene  image  at  the  seeker.  Our 
primary  concern  in  this  study  is  to  identify  potential  IR  obscurants  and  determine  their  extinction 
properties.  Our  initial  concern  is  the  3-5  pm  atmospheric  window;  however,  some  of  the  materials  also 
show  promising  characteristics  in  the  8-12  pm  window.  We  are  also  sensitive  to  potential  environmental 
and  health  related  issues,  ease  of  procurement  and  use,  and  availability  and  cost.  As  a  result,  all  initial 
candidates  are  readily  available  powders  with  acceptable  toxicological  and  optical  properties  that  are 
either  known  or  can  be  readily  determined  from  infrared  absorption  spectra. 

We  have  designed  and  constructed  an  experimental  apparatus  to  measure  the  extinction  spectra  of 
suspended  particles.  The  apparatus  consists  of  a  venturi  nozzle  to  lift  the  particles  and  mix  them  with 
nitrogen  carrier  gas.  The  gas/particle  stream  flows  up  a  flow  tube  equipped  with  optical  ports  for  Fourier 
Transform  Infrared  (FTIR)  interrogation.  Particle  mass  loading  is  calculated  from  the  measured  volume 
flow  rate  of  the  carrier  gas  and  the  mass  consumption  of  the  test  powder.  Using  this  apparatus,  we  have 
measured  the  mass  extinction  coefficients  for  brass  flakes,  325  mesh  boron  nitride  (BN)  powder  and  325 
mesh  boric  acid.  The  mass  extinction  coefficients  have  been  converted  to  volume  extinction  coefficients 
using  published  and  measured  packing  densities  of  the  individual  materials. 

In  parallel  with  the  experimental  efforts,  we  have  developed  a  series  of  codes  to  calculate  the 
absorbance  of  suspended  powders.  The  Mie  theory  provides  detailed  information  about  the  scattering  and 
absorption  processes  associated  with  individual  particles  in  the  obscurant  cloud.  The  quantities  calculated 
from  Mie  theory  (cross  sections,  phase  functions,  etc.)  can  be  related  to  experimental  observations  through 
the  radiative  transport  equation.  An  analytic  model  has  been  developed  to  describe  the  IR  absorbance  in  a 
cold  slab  geometry  of  suspended  particles. 

The  mass  extinction  coefficient  for  brass  powders  has  been  determined  from  IR  absorbance 
measurements  and  calculated  by  Mie  theory  using  experimental  optical  constants.  Our  experimental  and 
theoretical  results  for  the  mass  extinction  coefficient  agree  well,  but  are  lower  than  that  reported  in  the 
obscurant  literature  by  a  factor  of  15. 

We  performed  similar  experimental  and  theoretical  studies  on  BN  powder,  which  was  one  of  our 
initial  candidates.  We  found  that  the  spectra  in  the  literature  are  incorrect.  However,  we  have  derived 
good  optical  constants  to  model  its  absorption  and  scattering.  Although  BN  itself  does  not  have  strong 
resonances  in  the  3-5  pm  region,  it  does  exhibit  a  strong  resonance  at  7.4  \im  (OexF  0.27  m^/gm)  that  tails 
into  the  10  ^m  region.  This  may,  therefore,  be  useful  in  the  long  wave  IR  band. 
These  results  also  confirm  that  our  model  can  provide  useful  guidelines  for  predicting  accurate  cross 
sections  for  both  scattering  and  absorption  features  through  the  entire  IR,  provided  the  complex  index  of 
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refraction  and  particle  size  distribution  are  available.  We  are  confident  that  we  can  now  use  the  theory  to 
prefilter  candidate  materials,  and  provide  specific  guidelines  concerning  the  optimum  particle  sizes. 


To  move  the  absorption  into  the  3-5  p.m  region,  boric  acid  was  chosen  for  investigation.  The 
presence  of  an  OH  chemical  moiety  results  in  a  large  absorption  in  this  region,  and  boric  acid  can  be 
obtained  commercially  as  a  fine  powder  (325  mesh).  Although  the  spectra  appear  promising,  the  cross 
sections  are  too  low  (typically  -0.008  m^/gm)  for  it  to  be  an  effective  obscurant.  However,  it  could  be 
used  for  demonstration  purposes  or  proof  of  principle  tests.  We  are  optimistic  that  this  is  not  the  best 
system  available  in  this  region,  and  are  looking  at  several  alternatives. 

Extinction  properties  of  each  of  the  materials  tested  in  both  the  3-5  |xm  and  8-12  |Jm  bands  are 
shown  in  the  table  below.  The  fractional  extinction  shown  in  the  table  is  calculated  from  the  measured 
mass  extinction  coefficients  vs  wavelength,  the  mass  loading  calculated  from  the  packing  density  and  the 
assumed  charge  size  (70  cm^),  and  a  path  length  of  10  m.  The  resulting  wavelength-dependent  extinction  is 
multiplied  by  the  blackbody  curve  for  a  350K  source  in  the  3-5  pm  and  8-12  pm  regions.  The  values 
tabulated  are  the  average  extinction  in  these  two  windows. 


Fractional  Extinction  of  Radiation  from  a  350K  Blackbody  Resulting  from  70  cm^  of 
the  Materials  Investigated  in  This  Study  Dispersed  in  a  10-m  Diameter  Cloud 


Material 

Packing  Density 
(g/cm^) 

Band 

3-5  pm 

8-12  pm 

Brass  flakes 

2.6 

0.22 

0.22 

Boron  nitride 

0.666 

0.07 

0.05 

Boric  acid 

0.667 

0.02 

0.01 

The  methodology  and  approach  we  describe  is  both  sound  and  practical  and  can  be  applied  to  a 
number  of  problems.  We  propose  to  extend  this  technique  to  other  regions  of  the  spectrum,  including  the 
near-IR,  visible,  and  ultraviolet.  These  extensions  are  straightforward  applications  of  the  current 
experimental  setup,  with  some  additional  optics  to  cover  the  expanded  portions  of  the  spectrum. 
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OBSCURANTS  FOR  INFRARED  COUNTERMEASURES 


INTRODUCTION 

The  advent  of  imaging  seekers  in  the  3-5  gm  and  8-12  gm  bands  poses  a  new  challenge  in  the  area  of 
infrared  countermeasures.  Conventional  seekers  rely  on  the  power  received  in  a  particular  spectral  region 
or  regions,  but  they  are  not  sensitive  to  the  spatial  extent,  aspect  ratio,  or  velocity  of  the  perceived  target 
The  next  generations  of  seekers  will  be  based  on  arrays  that  will  provide  the  designer  with  a  wide  range  of 
options  for  discriminating  between  enemy  targets,  friendly  craft,  and  decoys.  The  simplest  discriminants 
involve  those  relating  to  spatial  extent  and  aspect  ratio.  Area-extensive  infrared  countermeasures  will  be 
necessary  to  defeat  these  new  seekers.  Based  on  work  by  the  Army,  one  effective  component  of  these  area- 
extensive  IR  countermeasmes  could  be  obscurants.  Our  primary  concern  in  this  study  is  to  screen  potential 
infrared  obscurants  as  to  their  extinction  properties,  ease  of  deployment,  chemical  inertness,  availability, 
cost,  and  environmental  effects. 

An  airborne  obscurant  modifies  electromagnetic  wave  propagation  between  a  potential  target  and  a 
posed  threat  or  observer.  The  radiation  from  the  target  can  be  removed  from  the  threat's  field  of  view  by 
scattering  and  absorption  processes.  This  reduction  in  radiance  is  classified  as  a  direct  transmission  loss. 
The  direct  transmission  loss  depends  on  the  obscurant's  mass  concentration,  spatial  distribution,  particle 
size  distribution  ,  particle  geometry,  as  well  as  its  bulk  optical  properties,  i.e.,  the  complex  index  of 
refraction.  The  spectral  signature  of  the  obscurant  cloud  is  affected  by  both  the  potential  target  and  the 
environment.  Forward  scattering  of  radiation  from  the  target  contributes  to  die  radiance  of  the  obscurant 
cloud.  Backward  scattering,  on  the  other  hand,  reduces  the  radiance.  Diffuse  multiple  scattering  within 
the  obscurant  cloud  can  introduce  radiation  from  point  sources  such  as  the  stm  or  extended  sources  such  as 
the  sky  or  natural  cloud  cover.  These  background  sources  contribute  to  the  obscurant  cloud's  radiance. 
Our  objective  in  this  program  is  to  find  materials  that  will,  with  some  combination  of  scattering  and 
absorption,  remove  spectral  radiance  in  specific  infiared  windows.  Initially,  our  concern  is  the  3-5  gm 
window;  however,  some  of  the  initial  candidates  also  show  promise  in  the  8-12  gm  window.  We  are  also 
sensitive  to  potential  environmental  and  health  related  issues,  ease  of  procurement  and  use,  and  availability 
and  cost.  As  a  result,  all  our  initial  candidates  will  be  readily  available  powders  with  acceptable 
toxicological  properties  and  optical  properties  which  are  either  known  or  can  be  readily  determined  from 
infrared  absorption  spectra. 

TECHNICAL  APPROACH 

This  project  encompasses  a  strong  interaction  of  experimental  and  theoretical  work.  The  major 
thrust  of  the  FY97  effort  was  the  identification  of  potential  materials  for  deployment  in  a  transmission  loss 
infrared  obscurant  mode.  The  materials  evaluated  included  both  metals  and  refractory  substances.  The 
metal  chosen  for  study  was  brass  in  the  form  of  flakes,  as  used  by  the  Army  in  their  obscurant  program. 
This  species  was  chosen  to  provide  a  baseline  for  both  the  experimental  and  theoretical  studies  to  follow. 
As  discussed  below,  at  this  point  we  do  not  consider  brass  a  viable  candidate  for  several 
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reasons.  However,  it  did  provide  us  with  the  opportunity  to  confirm  our  experimental  and  theoretical 
methodology.  The  other  materials  studied  included  boron  nitride  (BN)  and  boric  acid  (BA).  The  boron 
nitride  was  originally  chosen  because  of  an  apparent  absorption  peak  in  the  4  pm  region.  This  peak  was 
reported  by  Doll  (1994),  but  our  work  has  conclusively  shown  that  it  is,  in  fact,  not  due  to  BN.  This  issue 
will  be  addressed  in  the  discussion.  Chemically,  BN  is  a  refractory  material  that  is  relatively  inert.  It 
reacts  very  slowly  in  humid  air  to  produce  boron  oxides  and  ammonia.  Spectra  of  boric  acid  taken  by  us 
exhibit  relatively  strong  absorptions  in  both  the  3-5  pm  and  8-12  pm  regions.  Boric  acid  is  an  extremely 
weak  acid  (pH  =  5.2  for  a  O.IN  solution),  which  is  not  much  different  from  the  acidity  found  in  pure  water 
saturated  with  carbon  dioxide  (seltzer  water).  It  is  listed  as  an  irritant,  as  most  powders  are,  but  chemicall 
is  relatively  stable.  Neither  BN  nor  BA  are  considered  environmentally  hazardous. 

The  preliminary  theoretical  effort  was  directed  at  evaluating  spectral  scattering  and  absorption  cross 
sections  of  proposed  obscurants.  The  cross  sections  as  a  function  of  IR  wavelength  and  particle  size  were 
calculated  from  experimental  values  of  the  complex  index  of  refraction  using  the  Mie  theoiy  for  scattering 
from  spheres.  These  calculations  provided  information  on  the  relative  importance  of  scattering  and 
absorption  processes  for  a  proposed  material.  The  ratio  of  the  scattering  cross  section  to  the  sum  of  the 
scattering  and  absorption  (extinction)  cross  section  is  the  single-scattering  albedo.  If  the  albedo  is  small 
(«1  ),  the  spectral  transmittance  of  the  obscurant  cloud  can  be  readily  calculated  from  the  particle 
concentration  and  the  extinction  coefficient  using  the  exponential  Lambert-Beer  law. 

The  experimental  effort  measured  the  extinction  coefficients  of  proposed  obscurants,  A  gas  flow 
reactor  has  been  built  that  creates  a  suspension  of  particles  of  known  size  distribution  and  composition. 
The  extinction  of  IR  radiation  from  a  black-body  source  by  the  suspended  particles  is  measured  using  an 
Fourier  Transform  Infrared  (FTIR)  spectrometer.  Contributions  to  the  extinction  from  both  scattering  and 
absorption  are  measured,  and  the  theoretical  model  is  used  to  determine  the  albedo  and  particle  size 
distributions.  It  is  anticipated  that  the  particle  size  distributions  can  ultimately  be  tailored  to  obtain  th 
best  combination  of  scattering  vs  absorption  in  the  spectral  regions  of  interest.  In  this  way,  some  spectral 
tailoring  of  the  obscurant  can  be  accomplished.  In  all  cases,  spectral  measurements  are  compared  with  the 
theoretical  calculations  discussed  above.  With  guidance  from  these  calculations,  we  plan  to  iterate  on  the 
particle  size  distributions  and  chemical  nature  of  the  obscurant  to  find  the  optimum  material. 

EXPERIMENTAL  METHODS 

The  experiments  were  carried  out  in  an  aerosol  flow  tube  coupled  to  an  FTIR  spectrometer  (Fig. 
1).  Briefly,  powders  for  study  flow  up  an  aerosol  flow  tube  where  they  intersect  the  beam  of  the  FTIR 
spectrometer  (Mattson  Instruments,  Cygnus  100)  operating  with  an  external  detector.  Transmission 
spectra  are  recorded  with  and  without  particle  flow,  and  the  resulting  absorption  due  to  the  particles  is 
plotted  vs  wavelength  in  the  2-20  |Lim  region. 

Powders  to  be  tested  are  picked  up  using  a  venturi  nozzle  (Air-Vac  Engineering  Company,  Inc., 
model  TDl  lOH)  operating  with  compressed  nitrogen  as  the  feed  gas.  The  nozzle  is  designed  such  that  20 
standard  liters  per  minute  (SLPM)  of  feed  gas  produces  a  vacuum  flow  of  ~20  SLPM.  Powder  is  delivered 
to  the  vacuum  inlet  of  the  nozzle  by  a  rotating  platen  feeder  system.  It  consists  of  a  15-cm  diameter 
rotating  A1  platen  with  a  1  cm  wide  x  1.2  cm  deep  groove  at  a  radius  of  6  cm.  Powder  is  fed  into  this 
groove  opposite  the  venturi  feed  tube,  and  the  powder  level  is  adjusted  using  a  micrometer-controlled 
scraper.  The  platen  rotates  at  0.08  revolutions  per  second.  As  the  platen  rotates  and  the  powder 
approaches  the  venturi  nozzle,  it  is  lifted  by  the  vacuum  flow  through  a  6  mm  tube  into  the  nozzle  where  it 
is  mixed  with  the  feed  gas.  Most  experiments  are  conducted  with  a  total  flow  (feed  gas  plus  vacuum  flow) 
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of  40  SLPM.  Mass  loading  is  calculated  by  dividing  the  measured  mass  consumption  rate  by  this  flow 
rate.  Several  experiments  were  also  carried  out  at  flows  of  20  and  80  SLPM.  At  the  highest  flow,  window 
contamination  became  a  problem  and  these  data  were  not  used  in  the  analysis. 


to  hood 


Fig.  1  -  Schematic  diagram  of  the  apparatus  used  in  this  work 


After  exiting  the  venturi  nozzle,  the  combined  gas  and  particle  flow  enters  a  transition  funnel  that 
increases  the  area  of  the  flow  stream  and  lowers  the  linear  flow  velocity.  As  shown  in  Fig.  1,  this  funnel 
contains  the  main  flow  tube,  which  is  a  1 1 .4-cm  inside  diameter  (ID)  plexiglass  tube.  Two  sidearms  for  the 
optical  windows  are  located  on  opposite  sides  of  the  flow  tube,  35  cm  upstream  of  the  transition  funnel. 
The  windows  are  set  back  13  cm  from  the  wall  of  the  main  flow  tube.  A  nitrogen  purge  flow  of  1 .2  SLPM 
is  split  between  the  two  sidearms  to  ensure  that  the  powder  does  not  deposit  on  the  optical  windows.  After 
passing  through  the  probe  region,  the  particle  flow  is  exhausted  into  the  fume  hood. 

The  m  probe  beam  is  directed  out  of  the  spectrometer  using  a  spherical  mirror  witii  2-m  radius  of 
curvature.  This  results  in  a  nearly  collimated  beam  that  propagates  through  the  chamber  and  is  focused  on 
an  HgCdTe  detector.  The  spectra  are  recorded  in  single-beam  mode  from  400  to  5000  cm'’  (25  to  2  pm), 
with  a  spectral  resolution  of  1  cm"'.  Each  spectrum  recorded  is  the  summation  of  24  individual  scans. 
Particulate  extinction  spectra  are  calculated  by  normalizing  the  single  beam  spectrum  measured  during 
powder  flow  by  a  background  spectrum  collected  with  only  gas  flowing.  These  spectra  reflect  not  only 
particle  absorption  but  any  loss  due  to  scattering,  so  they  are  true  measures  of  the  extinction  coefficient  of 
the  particles. 
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EXPERIMENTAL  RESULTS 

Absorption  spectra  of  three  candidate  species  have  been  measured  over  the  2-14  pm  spectral  region 
under  several  different  flow  rates  and  particle  loadings.  Figure  2  is  a  representative  spectrum  of  brass.  All 
experimental  data  are  presented  in  absorbance  units  (absorbance  =  -log(I/Io).  Our  goal  in  measuring  the 
brass  extinction  was  to  have  a  system  to  use  as  a  benchmark  for  future  studies.  It  has  a  relatively 
unstructured  extinction  coefficient  throughout  the  IR,  possibly  indicating  the  dominance  of  scattering  over 
any  spectral  absorption  features.  As  such,  it  is  not  as  attractive  an  obscurant  as  some  of  the  other  species 
to  be  discussed.  It  has,  however,  been  used  by  the  Army.  When  this  program  started,  it  was  thought  to  be 
a  well  characterized  particulate  system  (Edwards  et  al.  1992).  However,  the  data  shown  are  not  consistent 
with  the  reported  mass  extinction  coefficients,  but  are  about  a  factor  of  15-20  lower.  Different  flow  rates 
and  particle  loadings  had  virtually  no  effect  on  the  mass  extinction  coefficients  determined  from  our 
experimental  data.  The  large  discrepancy  between  our  numbers  and  those  reported  by  Edwards  et  al.  led  us 
to  explore  some  of  the  potential  problems  areas  inherent  in  these  measurements,  i.e.,  multiple  scattering, 
particle  resonances,  and  scattering  anisotropy.  The  theoretical  model  and  results  are  discussed  in  the 
following  section,  however,  the  results  from  theoretical  calculations  support  our  measured  extinction 
coefficients.  At  this  time,  we  are  not  able  to  explain  the  Army’s  data  witii  any  certainty.  It  does  appear 
from  the  available  literature  that  the  combination  of  path  length  and  mass  loading  in  the  Army’s  experiment 
results  in  an  absorption  of  >98%  of  the  IR  radiation  incident  on  their  detector.  In  this  case,  it  is  quite 
possible  that  the  measurements  are  in  error  due  to  the  noise  floor  of  the  detector/amplifier  or  some 
saturation  effect.  Without  additional  details  on  the  exact  experimental  configuration,  it  is  impossible  to 
elaborate  further. 


Fig.  2  -  Experimental  absorption  spectrum  of  brass  flakes  at  135  g/m^  and  a  pathlength  of  1 1.4  cm 


The  prime  candidate  at  the  beginning  of  this  program  was  BN.  It  has  several  advantages  over 
brass.  It  is  readily  available  as  a  well-characterized  powder,  with  a  mass  weighted  distribution  peaking  at 
10  pm  diameter  particle  size.  It  is  non  toxic,  poses  no  environmental  hazard,  and  is  chemically  stable. 
According  to  Doll  (1994),  diere  is  a  strong  absorption  feature  in  the  3-4  pm  region.  Figure  3  shows  an 
absorption  spectrum  of  BN;  no  absorption  band  is  evident  in  this  spectral  region.  The  large  rising 
absorption  in  the  2  pm  region  is  thought  to  be  due  solely  to  scattering. 
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Fig.  3  -  Experimental  absorption  spectrum  of 325  mesh  BN  powder  at  98  g/m^  using  a  1 1.4  cm  path  length 

There  are  other  interesting  features,  however.  The  strong  absorption  feature  at  7.38  pm  carries 
well  into  the  8-12  pm  region  and  is  potentially  useful  against  seekers  in  this  region.  To  test  the  apparatus 
for  linearity,  a  series  of  experiments  at  various  flow  rates  and  BN  particle  loadings  was  undertaken.  These 
results  are  shown  in  Fig.  4  for  data  taken  in  three  spectral  regions.  It  is  evident  from  the  figure  that,  at  the 
flow  rates  being  used,  the  data  are  linear  over  at  least  an  order-of-magnitude  change  in  particle  loading. 
Thus  we  can  be  confident  that  these  data  can  be  extended  to  the  much  lower  concentrations  and  much 
longer  path  lengths  that  will  be  encountered  in  actual  deplojmient  tests. 


Fig.  4  ~  Measured  absorbance  of  BN  powder  at  three 
wavelengths  vs  mass  loading.  The  points  at  98  g/m^ 
correspond  to  volume  flow  rates  of  20  and  40  SLPM. 
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To  move  the  absorption  into  the  3-5  pm  region,  boric  acid  was  chosen.  The  presence  of  an  OH 
chemical  moiety  results  in  a  large  absorption  in  this  region,  and  boric  acid  can  be  obtained  commercially  as 
a  fine  powder  (325  mesh).  Figure  5  shows  a  spectrum  of  BA  particulates.  It  is  obvious  that  strong 
absorption  features  are  present  in  the  3-6  pm  and  7-1 1  pm  regions.  At  the  present  time,  we  do  not  have  an 
accurate  particle  size  determination  for  the  BA  samples  or  optical  constants  for  this  material,  but  these  can 
be  determined  and  optimized.  However,  more  work  must  be  done  to  determine  the  best  particle  size 
distribution  for  a  particular  spectral  region.  These  preliminary  data  are  very  encouraging  and  warrant 
additional  study  of  BA  as  an  obscurant. 


Fig.  5  -  Experimental  absorption  spectrum  of  boric  acid  powder  at  750  g/m^  and  1 1 .4  cm  pathlength 


THEORY 

Particle  Entrainment 

The  nature  of  the  gas-particle  flow  in  the  aerosol  flow  tube  described  above  depends  on  the  particle 
size  distribution,  the  particle  geometry  and  density,  the  gas  viscosity,  the  inlet  velocity,  and  the  chamber 
geometiy.  Particles  of  different  sizes  or  weight  can  be  sorted  by  the  flowing  gas.  Viscous  effects  of  the 
gas  are  important  near  the  tube  walls  and  at  the  particle  surfaces.  A  cylindrically  symmetric  stream 
function  model  (Zucrow  and  Hoffinan  1976)  has  been  developed  to  describe  the  gas  flow  within  the 
chamber.  The  stream  flmction  satisfies  the  continuity  equation  for  the  conservation  of  total  mass  in 
cylindrical  coordinates: 


V^'F  =  0, 
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where  p  is  the  gas  density  and  Vr  and  denote  the  radial  and  axial  velocity,  respectively.  The  radial 
component  of  velocity  is  nonzero  in  the  conical  base  region  of  the  test  chamber.  For  subsonic  (Mach 
number  «1)  or  incompressible  flow,  the  stream  function  is  a  solution  of  Laplace’s  equation  subject  to 
inlet  and  outlet  boundary  conditions.  The  velocity  components  are  calculated  from  the  spatial  derivatives 
ofT.  Since  the  tube  diameter  is  large  (11.4  cm),  the  stream  function  model  should  provide  a  reasonable 
approximation  to  the  flow  field  throughout  most  of  the  tube  volume.  However,  gas  viscosity  is  neglected  in 
this  model,  and  therefore  it  fails  to  predict  the  zero-velocity  (no-slip)  condition  near  the  tube  walls. 
Corrections  for  the  effects  of  gas  viscosity  on  the  flow  could  be  included  by  using  a  boundary-layer  model 
near  the  chamber  surfaces.  The  deposition  of  particles  on  the  chamber  walls  is  due  to  the  loss  of  “lift”  that 
occurs  within  the  boundary  layer.  Since  we  are  primarily  interested  in  the  bulk  flow  away  from  surfaces, 
the  boundary-layer  model  is  not  required.  The  key  results  of  the  stream  function  model  calculations  are 
given  in  Table  1  as  functions  of  the  inlet  flow  rate.  These  results  were  obtained  from  numerical  solutions 
of  the  Laplace  equation  using  a  finite-element  code  (PDEase  1995). 


Table  1  -  Results  of  the  Stream  Function  Model  Calculations 


Flow  rate  (liter/min) 

Inlet  Speed  (cm/s) 

Vz  (cm/s) 

Mass  Flow  Rate  (gm/s) 

20 

1052 

2.63 

0.374 

40 

2105 

5.26 

0.749 

80 

4210 

10.52 

1.497 

Over  the  range  of  inlet  flow  rates  considered,  flie  axial  velocity  (14)  drops  by  a  factor  of  400  within 
2.0  cm  of  the  inlet  exit.  Upon  leaving  the  conical  base  region,  the  model  predicts  a  uniform  plug  flow  up 
the  tube.  This  plug  flow  provides  “lift”  for  the  particles.  The  force  acting  upward  on  a  single  particle  of 
diameter  d  is  given  by  (Zhou  1993) 


where  Cj  is  the  drag  coefficient,  Vg  denotes  the  gas  stream  speed,  and  Vp  is  the  particle  speed.  In  the 
Stokes  regime,  the  drag  coefficient  for  a  spherical  particle  is  given  by 


r-  _  24v  24 

'•'d -jz 

\yg-yp\d  Re 

where  v  is  the  kinematic  viscosity  (N2/O.I586  cmVs  @  300K)  of  the  gas  and  Re  is  the  Reynolds  number. 
This  expression  for  the  drag  coefficient  is  valid  for  Reynolds  numbers  below  0.2  (creeping  flow).  The 
Reynolds  numbers  for  the  particles  considered  here  are  on  the  order  of  1  x  lO"^.  The  drag  coefficient  for  a 
disk  is  very  similar  to  that  of  the  sphere  (Chemical  Engineers'  Handbook  1973).  Thus,  the  analysis  is 
applicable  to  brass  flakes  modeled  as  falling  circular  disks. 

The  gravitational  force  (weight)  acting  downward  on  an  individual  spherical  particle  is  given  by 

^dn  ~  P  p§  g  ’ 
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where  pp  is  the  bulk  particle  density  (8.5  gm/cm^  for  brass)  and  g  is  the  gravitational  constant.  The 
magnitudes  of  these  forces  are  plotted  in  Fig.  6  as  a  function  of  brass  particle  diameter  for  various  flow 
rates.  The  lifting  force  is  a  linear  function  of  particle  diameter.  The  weight  increases  as  the  cube  of  the 
diameter.  The  forces  are  balanced  where  the  lines  intersect  the  weight  curve.  Particles  whose  diameters  lie 
to  the  left  of  the  intersections  are  accelerated  upward  in  the  gas  stream.  Particles  whose  diameters  are  to 
the  right  of  the  intersections  fall  downward.  The  absolute  velocity  of  a  particular  particle  depends  on  its 
terminal  velocity  as  well  as  the  gas  speed.  If  the  downward  terminal  velocity  is  equal  to  the  upward  gas 
speed,  the  particle  is  suspended,  having  an  absolute  velocity  of  zero.  These  calculations  assume  an 
equilibrium  flow  in  which  the  particle’s  diffiision  relaxation  time  is  very  short  compared  with  the  residence 
time  in  the  chamber.  The  chamber  residence  time  is  on  the  order  of  17  s  at  40  SLPM  and  the  diffusional 
relaxation  times  are  on  the  order  of  lO"^  to  10"^  s.  This  relaxation  time  implies  that  the  particles  respond 
quickly  to  the  flow  field  and  are  able  to  track  the  streamlines. 


brass  particle  diameter  (irm) 

Fig.  6  -  The  calculated  weight  of  spherical  brass  particles  and 
the  lift  force  they  experience  at  three  flow  rates  appropriate  to 
this  study 


The  calculations  reveal  that  at  40  SLPM  the  largest  spherical  brass  particle  we  can  “pick  up”  is 
about  14  pm  in  diameter.  If  the  brass  flakes  are  modeled  as  circular  disks  of  1.7  pm  diameter  with  a 
thickness  of  0.08  to  0.32  pm,  the  weight  is  on  the  order  of  1.0  x  10'®  dynes.  This  weight  is  four  orders  of 
magnitude  smaller  than  the  estimated  lift  forces  and  therefore,  all  flakes  should  be  entrained  in  the  main  gas 
stream. 

Radiative  Transfer  Models 

An  airborne  obscurant  material  modifies  electromagnetic  wave  propagation  between  a  potential 
target  and  a  posed  threat  or  observer.  The  radiation  from  the  target  can  be  removed  from  the  threat’s  field 
of  view  by  scattering  and  absorption  processes.  Each  of  these  processes  has  a  spectral  cross  section  (cm^) 
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associated  with  it,  which  is  a  function  of  particle  geomeby,  wavelength,  and  the  particle’s  complex  index  of 
refraction.  The  sum  of  the  spectral  absorption  cross  section  (Oa*)  and  the  spectral  scattering  cross  section 
(Osc)  is  the  spectral  extinction  cross  section  ((!«).  The  albedo,  denoted  by  ®,  is  the  ratio  of  the  scattering 
cross  section  to  tiie  extinction  cross  section 


(0= 


ab 


The  numerical  value  of  the  single-particle  albedo  indicates  the  importance  of  scattering.  An  albedo  of  one 
implies  that  the  obscurant  acts  purely  as  a  scattering  medium  in  which  no  absorption  or  emission  occurs  at 
the  given  wavelength.  An  albedo  of  zero  indicates  a  purely  absorbing  and/or  emitting  medium  with  no 
scattering.  Botti  the  albedo  and  cross  sections  depend  implicitly  on  the  radiation  wavelength  X  and 
polarization.  In  general,  both  absorption  and  scattering  processes  occur  simultaneously  within  the  cloud, 
and  the  albedo  lies  somewhere  between  zero  and  one.  The  loss  in  radiant  intensity  (I)  along  a  line  of  sight 
(LOS)  through  the  obscurant  cloud  is  classified  as  a  direct  transmission  loss  and  is  governed  by  the 
radiative  transfer  equation.  The  direct  transmission  loss  depends  on  the  obscurant’s  mass  concentration, 
spatial  distribution,  particle  size  distribution,  particle  geometry,  and  the  bulk  optical  properties,  i.e.,  th 
particle’s  complex  index  of  refraction. 

The  spectral  signature  of  the  obscurant  cloud  is  affected  by  both  the  potential  target  and  the 
environment.  Forward  scattering  of  radiation  from  the  target  contributes  to  tiie  obscurant  cloud’s  radiance 
(intensity).  Diffuse  multiple  scattering  within  the  cloud  can  introduce  radiation  from  point  sources  such  as 
the  sun  or  extended  sources  such  as  the  sky  or  natural  cloud  clover.  These  background  sources  all 
contribute  to  the  obscurant  cloud’s  radiance. 


The  Mie  theory  provides  detailed  information  about  the  scattering  and  absorption  processes 
associated  with  individual  particles  in  the  obscurant  cloud.  The  quantities  calculated  from  Mie  theory 
(cross  sections,  phase  functions,  etc.)  can  be  related  to  experimental  observations  through  the  radiative 
transport  equation.  An  analytic  model  is  now  developed  to  describe  the  ER.  absorbance  in  a  cold  slab 
geometry  of  suspended  particles  (Fig.  7). 


Fig.  7  -  Schematic  illustrating  the  slab  geometry  used  in  the  absorbance  calculations 


The  radiative  transfer  equation  for  an  isotropic  slab  of  thickness  L  with  azimuthal  symmetty  is 
given  by  (Ozisik  1972) 
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^  -1 


where  )J.  is  the  cosine  of  the  angle  0  between  the  direction  of  the  radiation  intensity  and  the  positive  OX 
axis  through  the  slab,  Np  is  the  number  of  particles  per  unit  volume  (1/cm^),  ft'  denotes  the  cosine  of  the 
in-scattering  angle  0%,  and  p(|X,|J.')  is  the  phase  function.  The  phase  function,  which  depends  on 

wavelength  X,  the  particle  diameter,  and  the  material’s  optical  constants,  provides  information  about  the 
directional  distribution  of  scattered  radiation.  For  the  slab  geometry  considered  here,  the  phase  function 
describes  scattering  from  the  p.'  direction  (a  cone  of  incident  directions  making  an  angle  of  0's  with  the  OX 
axis)  into  the  p  direction  (a  cone  of  scattering  directions  making  an  angle  of  0  with  the  OX  axis).  In  the 
laboratory,  the  absorbance  measurements  are  done  at  0  =  0  (p  =  1) ,  but  the  in-scattering  angle  varies  from 
0  to  180  deg.  If  the  particle  size  distribution  is  known  within  the  cloud,  the  extinction  coefficient  (y  = 
Np[Oab  +<7sc])  and  phase  function  are  calculated  by  a  suitable  average  over  the  given  particle  distribution. 

The  first  term  on  the  right-hand  side  of  the  radiative  transfer  equation  represents  the  loss  in  radiant 
intensity  I  due  to  absorption  and  scattering  out  of  the  line  of  sight  (LOS).  The  last  integral  term  represents 
the  gain  in  radiant  energy  by  the  beam  along  the  LOS  due  to  in-scattering  of  radiation  by  the  surrounding 
obscurant.  The  in-scattering  term  decreases  the  obscurant  cloud’s  absorbance.  The  efficiency  of  in¬ 
scattering  to  increase  the  radiant  intensity  depends  on  the  phase  function  and  the  magnitude  of  the  albedo  ©. 
The  phase  function  is  often  represented  by  expansions  in  Legendre  polynomials  (Pn(x))  of  the  form 

71=0 

where  ao  =  1.  The  relation  between  the  slab  variables  (|l',|Lt)  and  the  single-particle  phase  function 
calculated  from  Mie  theory  can  be  obtained  from  solid  geometry. 

If  the  obscurant  particles  are  assumed  to  be  spherical,  homogeneous,  and  isotropic,  the  Mie 
scattering  theory  can  be  used  to  calculate  the  spectral  absorption  and  scattering  cross  sections,  the  albedo, 
and  the  phase  functions  for  the  specified  particle  size.  In  general,  the  obscurant  particles  are  not  all  the 
same  size  (monodispersed),  and  the  Mie  calculations  would  have  to  be  done  with  a  given  particle 
distribution  function.  Since  we  have  no  information  on  the  particle  size  distribution  for  the  brass  flakes,  Np 
denotes  the  number  of  monodispersed  particles  per  cubic  centimeter. 

When  the  albedo  is  less  than  one,  the  radiative  transfer  equation  can  be  solved  for  the  spectral 
intensity  by  a  perturbation  series  in  co 

I{x,X)=Iq  +co*Ii  +00^ -12  +... 


Substitution  of  this  series  solution  into  the  radiative  transfer  equation  yields  in  lowest  order  the  Lambert- 
Beer  law 
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dJ(,_  iVp(ga6  +  g^)-Jo  _-Y-Jo 
dx  ^  n 


kU 


|U,)  =  /(0)-exp 


-y-  X 


=-logio 


I{0) 


>1 

J 


where  1(0)  denotes  the  initial  spectral  intensity  and  A(X,  x,  [i)  is  the  Lambert-Beer  spectral  absorbance. 
The  product  Np(Gab  ^Osc)  =  Y  defines  the  spectral  extinction  coefficient  in  units  of  cm'\ 

When  the  phase  function  is  represented  by  a  simple  analytic  expression,  e.g.,  the  Legendre 
expansion,  the  solution  Iq  can  be  substituted  into  the  in-scattering  integral  and  the  first-order  correction  Ii 
obtained  in  closed  form.  The  simplest  phase  function  is  for  isotropic  scattering  (p  =  1)  for  which  the 
spatial  distribution  of  scattered  radiation  about  the  scattering  center  (particle)  is  uniform.  The  phase 
function  for  Rayleigh  scattering  of  unpolarized  radiation  is  p(0^)  =  3/4(1  +  cos^  9^) ,  where  0p  denotes 

the  particle’s  single-scattering  polar  angle.  9p  is  the  angle  between  the  in-scattering  direction  and  the  out- 
scattering  direction  in  the  particle’s  frame  of  reference.  The  Rayleigh  scattering  distribution  has  equal 
“lobes”  in  the  forward  and  backward  direction  about  the  scattering  center.  Figure  8  presents  calculated 
phase  functions  for  brass  spheres  of  different  diameters  at  =  4.0  \xm  and  10.0  pm.  At  =  4.0  pm,  the 
Mie  scattering  is  veiy  anfsotropfc  for  larger  particles.  Forward  scattering  (0  <  05c<  90)  is  favored  over 
backward  scattering  for  the  4.83  pm  and  2.90  pm  diameter  particles.  The  phase  function’s  angular 
dependence  is  weak  for  particle  diameters  <1.5  pm.  In  the  3-5  pm  band,  the  small  particle  phase  functions 
are  almost  isotropic.  At  ^  =  10.0  pm,  the  phase  functions  for  the  smaller  particles  are  very  similar,  while 
the  larger  particle  diameters  now  display  stronger  backscattering.  If  the  extinction  is  dominated  by 
scattering,  the  measured  absorbance  depends  strongly  on  the  magnitude  of  the  in-scattering  integral. 


Fig.  8  -  Calculated  phase  functions  of  brass  spheres  of  different  diameters 
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Higher-order  corrections  to  the  Lambert-Beer  Law  can  be  obtained  from  perturbation  theory,  but 
they  are  cumbersome  and  usually  require  numerical  evaluation  of  the  in-scattering  integrals.  Rather  than 
deal  with  the  full  phase  function  over  all  wavelengths,  it  is  sometimes  more  convenient  and  useful  to 
consider  some  moment  of  the  angular  distribution  that  describes  the  relative  forward-to-backward 
scattering  ratio.  The  asymmetiy  factor  g  is  defined  as 


g  = 


—  f p*cos9  d(cos0). 


The  value  of  granges  from  -1  (maximum  backward  scattering)  to  0  (isotropic  scattering)  to  +1  (maximum 
forward  scattering).  The  calculation  of  this  quantity  as  a  function  of  wavelength  provides  insight  into  the 
spectral  variation  of  the  phase  function.  Figure  9  gives  calculated  asymmetry  parameters  for  brass  spheres 
of  various  diameters  as  a  function  of  wavelength.  The  larger  diameter  particles  tend  to  scatter  the  3-5  pm 
radiation  in  the  forward  direction  (g  >  0)  and  the  8-12  pm  radiation  in  the  backward  direction  (g  <  0)  for 
diameters  <  2.9  pm.  The  smaller  diameter  particles  act  as  backscatterers  over  the  entire  spectral  region. 


Fig.  9  -  Calculated  asymmetry  parameters  for  brass  spheres  of 
various  diameters  as  a  function  of  wavelength 


Scattering  effects  should  be  incorporated  into  the  radiative  transfer  model  for  accurate  calculation 
of  absorbance  and  interpretation  of  experimental  measurements.  The  presence  of  the  in-scattering  integral 
term  in  the  radiative  transfer  equation  presents  the  primaiy  difficulty  in  solving  the  equation  for  radiant 
intensity.  The  way  in  which  tiiis  term  is  treated  distinguishes  the  various  methods  of  solution  (Brewster 
1992).  One  of  the  simplest  techniques  is  the  two-flux  method.  The  two-flux  method  groups  all  scattering 
into  either  the  plus  or  minus  x  direction  along  a  LOS.  Assuming  axial  symmetry,  the  equation  of  radiative 
transfer  along  the  LOS  is  replaced  by  two  linear  equations 
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dx 


-  J  /(G)  •  r(x,  0)  •  sin  GcfG 

^  0 


/ 


y=Np(c,h+a,c) 


where  I  and  I '  denote  the  spectral  radiation  intensity  at  the  distance  x  along  the  LOS  in  the  directions  of 
increasing  and  decreasing  x,  respectively,  I'{x,  G)  is  the  radiation  intensity  reaching  a  scattering  center  at 

the  angle  0  with  respect  to  the  LOS,  and  /(0)  is  the  scattering  angular  distribution  at  the  angle  0.  In  the 
single-scattering  approximation,  i(0)  is  the  phase  fimction  obtained  from  Mie  theory.  The  integral  in  the 
above  equations  is  now  replaced  by 

J  /(0)  •  I'ix,  0)  •  sin  0d0  =  It  (x)  ■  J  /(0)  •  sin  0d0  -t-  /;  (x)  ■  f  f(e)  ■  sin  0d0 

0  0  71/2 


71/2 

J  /(G)  •  /'(a,  0)  ♦  sin  6d6 

__o _ 

7:/2 

J  /(0)-sined0 

0 


7C 

f  /(0)sin  6d6 

n/2 


where  It(x)  and  (x)  correspond  to  the  forward  and  backward  scattering  intensity,  respectively.  The 

two-flux  model  assumes  that  It(x)=  I*(x)  and  I;{x)  =  I-{x).  This  assumption  reduces  the  system  of 
integro-differential  equations  to  coupled  ordinary  differential  equations 

= -lVp(c,j  +  c^)  •  r  +  Npa,,a  ■  I*  +  NpOj  ■  r  r{o)=  1(0) 
^=^pi<^ab+<^sc)-i~-  Npc^a  ■  r  -  iVp<j,,p  •  r  r  (l) = o 


ti/2  n 

where  a  =  27t  J  /(0)*sin0d0  and  P  =  27C  J  /(G)*  sin  0cf0.  Note  that  the  a  and  p  terms  have  opposite 

0  7C/2 

signs  from  the  loss  terms  in  each  equation  and  increase  the  radiation  intensity  along  the  LOS.  The  two 
equations  on  the  right  provide  initial  conditions  for  the  system’s  solution.  The  condition  /  (0)  =  1(0) 
represents  the  incident  radiation  from  the  IR  source.  The  second  condition  I  ~  (L)  =  0  represents  a 
transparent  boundary  on  the  detector’s  side  of  the  slab.  These  equations  are  similar  to  the  classical 
Kubelka-Munk  equations  (Klein  1964)  for  diffuse  radiation. 

Both  a  and  P  are  functions  of  the  radiation  wavelength  X.  To  solve  the  equations  at  all  desired 
wavelengths  of  interest,  the  Mie  theory  must  be  used  to  evaluate  both  the  cross  sections  and  phase 
functions  at  each  wavelength.  A  FORTRAN  subroutine  has  been  developed  at  NRL  to  compute  the 
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integrals  a  and  P  at  each  wavelength  from  the  Mie  code  output.  The  solution  of  this  coupled  linear  system 
of  equations  gives  an  expression  for  the  spectral  absorbance  of  the  slab  geometry 

.  ,  ,  r 

‘ogio  -  *°g'Orcosh(rL)+i),sinh(rL)’ 

where  t{  =  Np{G^i,+G^{l-a))md  r  =  -Jb^  +  NpGlJ^^  .  If  a  and  p  are  zero,  then  r=  bi  and  this 

expression  for  the  spectral  absorbance  reduces  to  the  Lambert-Beer  law.  Note  that  the  quantities  F  and  bi 
are  dependent  on  wavelength  X  through  the  integrals  a  and  p. 

Brass  Extinction  Cross  Section 

The  Mie  theory  calculation  of  the  extinction  cross  section  requires  the  particle  geometry  and 
component  material’s  optical  constants  as  input  data.  Optical  constants  for  brass  were  taken  from 
Querry’s  experimental  work  (Querry  1985).  Figure  10  shows  the  complex  index  of  refraction  for  brass  as 
a  function  of  wavelength.  The  brass  used  in  this  experiment  came  from  an  M76 IR  screening  grenade.  The 
brass  powder  for  these  grenades  is  provided  by  Wolstenholme  Bronze  Powders  (WBP)  in  two  grades.  The 
differences  between  the  two  grades  are  the  percentage  of  stearic  acid  added  as  a  milling  aid  and  die  particle 
size  distribution.  The  WBP  high  grease  MD7000  brass  powder  has  a  mean  diameter  of  4.17  pm,  with  90% 
of  the  powder  less  than  7.25  pm.  The  Rich  MD7000  powder  has  a  mean  diameter  of  4.83  pm,  with  90% 
of  the  powder  less  than  10.26  pm  (Edwards  et  al.  1992). 


Fig.  10  -  The  real  and  imaginary  parts  of  the  index  of 
refraction  for  brass  in  the  3  -  13  pm  region 


The  extinction  cross  section,  the  scattering  cross  section,  and  the  absorption  cross  section  are 
calculated  as  a  function  of  wavelength  X  (3.0-14.0  pm)  using  an  NRL-modified  Mie  code  (Dave  1968). 
The  Mie  calculations  provide  the  single-particle  albedo  ©,  phase-functions,  and  asymmetry  parameters  g 
as  functions  of  X.. 


Obscurants  for  Infrared  Countermeasures 


15 


Tables  2  and  3  present  some  results  of  the  Mie  computations  for  nonpolarized  radiation  in  the  3-5 
pm  and  8-12  pm  bands.  All  of  these  results  are  numerical  averages  of  the  Mie  calculations  for  the 
specified  particle  diameter  d  over  the  respective  spectral  window.  The  number  of  monodispersed  particles 
per  cubic  centimeter  Np  is  based  on  the  experimental  mass  loading  of  135  g/m^ . 


Table  2  -  Results  of  Mie  Calculations  for  Brass  in  the  3-5  pm  Band 


d(pm) 

Np  (cm'^) 

Oext  (m’/g) 

CTjc  (cm’) 

cra6(cm’) 

CO 

g 

0.50 

2.43x10* 

0.061 

2.47x1  O''* 

9.03x10’" 

0.6396 

-0.2755 

0.90 

4.16x10’ 

0.471 

6.97x10'’ 

4.44x10’” 

0.9198 

-0.2253 

0.954 

3.54x10’ 

0.245 

8.95x10’’ 

5.14x10’” 

0.9307 

-0.2140 

1.07 

2.50x10’ 

0.269 

1.39x10-* 

6.69x10’'° 

0.9461 

-0.1873 

1.50 

8.99x10* 

0.267 

3.88x10-* 

1.33x10’’ 

0.9667 

-0.0393 

2.30 

2.49x10* 

0.177 

9.31x10-* 

2.71x10’’ 

0.9717 

0.2211 

2.90 

1.24x10* 

0.140 

1.48x10-’ 

1.48x10’’ 

0.9729 

0.3345 

3.70 

5.99x10* 

0.109 

2.39x10’’ 

6.39x10’’ 

0.9742 

0.3952 

4.83 

2.69x10* 

0.082 

4.02x10’’ 

1.02x10-* 

0.9750 

0.4436 

Since  the  albedo  ©  in  Table  2  is  close  to  one  for  most  of  the  particle  sizes  considered,  the  observed 
spectral  extinction  in  the  3-5  pm  band  is  primarily  due  to  scattering.  Particles  whose  diameters  are  below 
1 .5  pm  tend  to  backscatter  the  IR  radiation  (g  <  0).  Particles  with  diameters  on  the  order  of  1 .5  pm  tend  to 
be  isotropic  scatters  (g  =  0),  while  particles  with  diameters  above  1.5  pm  tend  to  scatter  the  IR  radiation  in 
the  forward  direction.  These  results  are  consistent  with  the  phase  fimctions  plotted  in  Fig.  8. 


Table  3  -  Results  of  Mie  Calculations  for  Brass  in  the  8-12  pm  Band 


d(pm) 

Oe;ff(m’/g) 

OscCcm’) 

<7a6(cm’) 

0) 

8 

0.50 

0.011 

4.80x10’'’ 

5.46x10’" 

0.0750 

-0.3036 

0.90 

0.012 

1.75x10’” 

2.06x10’'° 

0.4112 

-0.3283 

0.954 

0.013 

2.5010’” 

2.3510’” 

0.3719 

-0.3282 

1.07 

0.015 

4.95x10’'° 

3.94x10’'° 

0.5617 

-0.3268 

1.5 

0.031 

3.96x10’’ 

6.87x10’” 

0.8108 

-0.3096 

2.30 

0.084 

4.33x10'* 

2.14x10’’ 

0.9441 

-0.2500 

2.90 

0.108 

1.14x10’’ 

3.70x10’’ 

0.9665 

-0.1888 

3.70 

0.106 

2.33x10’’ 

5.85x10’’ 

0.9754 

-0.0780 

4.83 

0.083 

4.08x10’’ 

9.09x10’’ 

0.9782 

0.1050 

The  average  mass  extinction  coefficients  Gext  are  comparable  for  the  0.50  pm  and  4.83  pm 
diameter  particles.  However,  the  extinction  process  for  the  smaller  particle  involves  significant  absorption 
(©  =  0.64).  These  calculations  demonstrate  the  difficulty  of  obtaining  unique  particle  sizes  from  extinction 
measurements.  The  extinction  coefficient  y  for  the  0.5  pm  particle  is  0.082  cm''.  The  extinction  coefficient 
for  the  4.83  pm  diameter  particle  is  0.108  cm''.  In  the  small  particle  limit,  it  is  possible  to  obtain  a  particle 
diameter  that  would  give  an  extinction  coefficient  identical  to  that  of  the  larger  particle.  The  calculated 
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mass  extinction  coefficients  plotted  in  Fig.  11  as  a  fimction  of  particle  diameter  readily  demonstrate  this 
ambiguity.  Any  horizontal  line  intersecting  the  y  axis  defines  a  imique  extinction  coefficient,  but  the  line 
intersects  the  calculated  extinction  curve  at  two  points,  corresponding  to  a  small  particle  diameter  and  a 
large  particle  diameter.  The  solution  is  not  unique.  In  the  large  and  small  particle  limits,  the  extinction 
coefficients,  albedo,  and  phase  function  become  independent  of  the  particle  size  distribution. 


particle  diameter  (pm) 
(a)  3-5  pm  band 


particle  diameter  (pm) 
(b)  8-12  pm  band 


Fig.  1 1  —  Calculated  mass  extinction  coefficients  of  brass  particles  as  a  function  of  particle  size. 


Table  3  shows  that  scattering  is  dominant  for  particles  larger  than  1.5  pm  in  the  8-12  pm  band  and 
absorption  is  significant  for  particles  below  0.50  pm.  The  mass  extinction  coefficient  in  this  band  is 
generally  smaller  than  that  of  the  3-5  pm  band.  The  extinction  coefficient  of  the  4.83  pm  diameter  particle 
is  almost  identical  in  both  bands,  but  the  asymmetry  parameters  are  quite  different.  In  the  3-5  pm  band, 
the  4.83  pm  particle  tends  to  scatter  radiation  in  the  forward  direction,  but  the  scattering  is  roughly 
isotropic  in  the  8-12  pm  band. 


Brass  Absorbance  Calculations 

The  spectral  absorbance  can  be  calculated  from  the  computed  extinction  cross  section,  the  particle 
loading  and  the  geometric  path  length  of  the  chamber  (L  =  11.4  cm)  using  either  the  Lambert-Beer  law 
model  or  the  two-flux  model  described  above.  These  calculations  depend  on  the  assumed  particle  size 
distribution.  Since  the  relative  size  distribution  is  unknown,  it  is  assumed  to  be  monodisperse  for  each 
calculation.  Figure  12  shows  the  calculated  absorbance  based  on  the  Lambert-Beer  law  model  for  various 
spherical  diameters  of  brass  particles.  The  experimental  mass  loading  is  used  to  calculate  iVp  as  a  function 
of  particle  diameter  (Table  2).  The  calculated  absorbances  decrease  with  increasing  wavelength  for 
particle  diameters  below  3.7  pm.  This  result  is  consistent  with  the  calculated  extinction  coefficients  given 
in  Fig.  11.  At  a  particle  diameter  of  4.83  pm,  the  extinction  cross  section  in  the  3-5  pm  band  is 
independent  of  wavelength.  In  the  8-12  pm  band,  the  extinction  cross  sections  are  slightly  lower  and 
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roughly  equal  to  each  other.  The  experimental  absorbance  A  is  quite  flat  over  the  measured  spectral  range 
(A  =  0.46).  The  cross  section  calculated  for  the  4.83  pm  diameter  particle  gives  a  flat  absorbance.  The 
small  diameter  particle  solution  does  not  give  a  flat  absorbance  profile. 


wavelength  (gm) 

Fig.  12  -  Calculated  absorbance  for  a  loading  of  135  g/m^ 
and  pathlength  of  11.4  cm  for  various  spherical  diameters 
of  brass  particles  based  upon  the  Lambert-Beer  law  model 

The  Mie  calculations  indicate  that  the  extinction  by  brass  particles  is  dominated  by  scattering 
throughout  the  spectral  range  (3-12  pm).  The  absorbance  was  also  calculated  from  the  two-flux  model 
using  the  same  Mie  cross  sections.  The  Mie  code  provided  the  spectral  phase  functions  for  various  sphere 
diameters  at  each  wavelength.  These  phase  functions  were  averaged  over  the  scattering  angle  to  give  the 
spectral  parameters  ct  and  P  defined  above.  Figure  13  shows  the  two-flux  absorbance  calculations  for  the 
same  particle  diameters  as  Fig.  12.  The  inclusion  of  the  in-scattering  term  in  the  radiative  transport  model 
lowers  the  absorbance  by  30-50%.  Over  the  3-5  pm  band,  the  4.83  pm  diameter  particle  has  a  strong 
forward  scattering  component  (Fig.  8),  and  the  a  parameter  is  0.75  and  the  P  parameter  is  0.23.  The 
decrease  in  calculated  absorbance  over  this  band  from  0.54  for  the  Lambert-Beer  law  model  to  0.15  for  the 
two-flux  model  is  due  to  the  in-scattering  of  radiation  toward  the  detector.  In  the  8-12  pm  band,  a  =  p  = 
0.5  and  the  in-scattering  terms  cancel.  The  absorbance  approaches  the  Lambert-Beer  model. 

In  Fig.  14,  the  extinction  coefficients  derived  from  the  experimental  data  are  compared  to  those 
calculated  above.  Mass  extinction  coefficients  are  defined  as  =  2.303  •  AJL-C,  where  Ax  is  the 

measured  absorbance  at  wavelength  X,  L  is  the  path  length,  and  C  is  the  mass  loading  in  g/m^.  The  volume 
extinction  coefficients  are  calculated  by  multiplying  by  the  material  packing  density,  either  from  the 
literature  or  measured  in  this  work.  The  calculated  cross  section  from  the  4.83  pm  diameter  Mie  solution 
was  used  for  both  the  Lambert-Beer  model  and  the  two-flux  model.  The  particle  distribution  is  assumed  to 
be  monodispersed.  The  particle  diameter  of  4.83  pm  gives  an  almost  flat  Lambert-Beer  law  absorbance 
over  the  3-14  pm  wavelength  region,  which  lies  above  the  experimental  curve.  The  effects  of  the  in¬ 
scattering  term  on  the  absorbance  is  shown  in  the  same  figure  by  the  two-flux  model.  The  absorbance  is 
reduced  by  a  factor  of  3.7  in  the  3-5  pm  band.  This  result  is  readily  understood  from  the  asymmetry 
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Fig.  13  -  Calculated  absorbance  for  the  same  particle 
diameters  as  Fig.  12  using  the  two-flux  model 


parameter  g=  0.444  of  Table  2  and  the  phase  function  depicted  in  Fig.  8.  In  the  3-5  pm  band,  the  4.83  pm 
diameter  brass  powder  scatters  radiation  into  the  forward  direction.  In  the  8-12  pm  band,  the  asymmetry 
parameter  is  0.105  and  the  forward  scattering  component  of  the  intensity  is  less.  Thus,  the  reduction  in 
absorbance  is  not  as  great.  The  experimental  absorbance  lies  between  these  two  theoretical  absorbance 
models.  Only  single  scattering  processes  are  considered  in  the  two-flux  model.  The  inclusion  of  multiple 
scattering  would  increase  the  absorbance  in  the  3-5  pm  band.  The  effects  of  particle  size  distribution  on 
the  calculation  have  not  been  fully  explored. 


Fig.  14  -  Comparison  of  the  extinction  coefficients  calculated 
from  the  measured  absorbance  of  brass  flakes  with  the  results 
of  the  two  models  assuming  4.83  pm  diameter  particles 
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The  mass  extinction  coefficient  calculated  for  brass  particles  with  a  mean  diameter  of  4.83  |xm  in 
both  the  3-5  |im  band  and  the  8-12  |im  band  is  0.082  m^/g.  The  values  reported  by  Edwards  et  al,  are  1.13 
m^/g  (3-5  |Lim)  and  1.34  m^/g  (8-12  jim).  These  latter  values  are  higher  than  ours  by  a  factor  of  14-16. 
However,  our  experimental  results  are  consistent  with  the  theoretical  Mie  calculations  and  the  mean 
particle  size  provided  by  WBP.  The  agreement  of  the  calculated  brass  absorbance  and  the  experimental 
absorbance  provides  strong  evidence  that  our  mass  extinction  coefficient  is  correct. 

Boron  Nitride  Calculations 

Boron  nitride  has  four  crystal  structures  (Doll  1994):  hexagonal,  cubic,  wurtzitic,  and 
rhombohedral.  In  the  hexagonal  form,  it  is  a  slippery  white  solid  with  a  layered  structure  very  similar  to 
graphite  (Cotton  and  Wilkinson  1972).  At  high  temperature  and  pressure,  the  hexagonal  form  is  converted 
to  a  cubic  form  with  a  diamond  structure.  The  cubic  form  is  very  hard  and  will  scratch  diamond.  BN  is 
stable  in  dry  air,  but  is  slowly  hydrolyzed  by  water  to  produce  boron  oxides.  The  hexagonal  form  is 
commercially  available  from  The  Carborundum  Company,  Amherst,  NY.  Several  references  on  the 
dielectric  properties  of  the  hexagonal  form  are  available  in  the  open  literature  (Geick,  et  al.  1966;  Hoffman 
et  al.  1984).  The  optical  constants  of  cubic  BN  over  a  photon  energy  range  of  2.5-20  eV  has  been  reported 
by  Miyata  et  al.  (1989).  Geick  et  al.  (1966)  have  performed  a  Kramers-Kronig  analysis  of  the  IR  reflection 
spectra  of  hexagonal  BN  from  100-3000  cm“^  The  analysis  provides  the  complex  dielectric  fimctions  £(X) 
for  IR  radiation  polarized  parallel  and  perpendicular  to  the  hexagonal  c  axis.  These  dielectric  functions  are 
represented  by  a  classical  oscillator  fit  with  resonances  at  12.8  pm  (E  parallel  to  c)  and  12  pm  (£ 
perpendicular  to  c)  over  a  spectral  range  5  -33  pm.  The  absorption  spectrum  exhibits  several  minor  peaks 
that  are  probably  due  to  2-phonon  combination  processes.  Hoffman  et  al.  (1984)  have  performed  optical 
reflectance  studies  of  the  c  face  of  BN  over  a  photon  energy  range  of  0.045-10  eV  and  obtained  dielectric 
functions  that  are  quite  similar  to  those  of  Geick  et  al.  However,  we  question  the  accuracy  of  their  data 
over  the  3-5  pm  range.  They  report  a  large  absorption  feature  that  we  do  not  observe  in  IR  transmission 
and  cannot  explain  on  theoretical  grounds.  Their  sample  may  have  been  contaminated. 

The  complex  index  of  refraction  m  =  n  - 1.  •  k  can  be  calculated  from  the  square  root  of  8(X),  The 

optical  constants  of  hexagonal  BN  have  been  calculated  over  the  5-12  pm  range  using  the  dielectric 
functions  provided  by  the  Geick  et  al.  oscillator  models.  Because  these  data  are  not  published,  it  is 
tabulated  in  Table  4  for  reference. 


The  Mie  theory  is  applicable  to  a  homogeneous  spherical  particle.  If  the  particles  are  small 
compared  to  the  wavelength,  the  exact  Mie  theory  of  scattering  reduces  to  an  electrostatic  approximation 
(Rayleigh  limit).  If  the  particle  is  a  small,  homogeneous,  isotropic  sphere,  the  incident  radiation  presents  a 
uniform  electric  field  to  the  particle  and  induces  an  electric  polarization  that  results  in  a  dipole  parallel  t 
the  incident  field.  Since  the  sphere’s  dielectric  function  is  a  scalar  quantity,  its  polarizability  does  not 
depend  on  its  orientation.  In  the  Rayleigh  limit,  the  cross  sections  for  extinction  and  scattering  are  directl 
related  to  the  polarizability.  Thus,  for  a  spherical  particle  composed  of  isotropic  material,  the  cross 
sections  are  independent  of  the  sphere’s  orientation.  Hexagonal  BN  is  an  anisotropic  material.  The 
polariziability  of  hexagonal  BN  depends  on  the  orientation  of  the  incident  field  relative  to  the  c-axis  (Table 
4).  The  dipole  induced  within  a  spherical  particle  of  BN  by  the  incident  field  is  not,  in  general,  parallel  to 
that  field.  The  polarizability  is  a  tensor.  Determination  of  the  extinction  and  scattering  cross  sections 
under  these  conditions  requires  some  care.  The  polarizability  of  an  anisotropic  sphere  can  be  calculated 
exactly  in  the  electrostatic  approximation.  The  average  extinction  cross  section  of  a  collection  of  randomly 
orientated,  sufficiently  small,  anisotropic  spheres  is  (Bohren  and  Huffman  1983) 
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r  j  r 

^  ^ _ ^parallel  ,  ^  ^perpendicular 

m  3  3  . 

where  C,  denotes  the  cross  section  of  an  isotropic  sphere  with  i-th  polarization.  This  formula  is  valid  for 
particles  that  are  small  compared  to  X.  In  hexagonal  BN,  the  two  principal  directions  perpendicular  to  the 
c  axis  are  equivalent  (2/3).  Fig.  15  compares  the  average  mass  extinction  coefficient,  <Cea>,  with  the 
experimentally  determined  value.  The  mass  extinction  coefficients  were  calculated  at  specific  wavelengths 
using  the  derived  optical  constants  (Table  4)  and  a  particle  size  distribution  function  provided  by  The 
Carborundum  Company. 


Table  4  -  Calculated  BN  Complex  Index  of  Refraction 


Electric  field  perpendicular  to  caxis 

Electric  field  parallel  to  c  axis 

n 

k 

n 

k 

6.0 

1.037 

0.099 

1.426 

0.183 

6.5 

0.215 

1.454 

1.671 

1.690 

7.0 

0.603 

3.829 

2.740 

0.309 

7.5 

6.212 

1.113 

2.410 

0.088 

8.0 

4.012 

0.167 

2.216 

0.044 

8.5 

3.459 

0.074 

2.197 

0.029 

9.0 

3.202 

0.045 

2.136 
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Fig.  15  -  Comparison  of  experimentally  determined  mass 
extinction  coefficients  for  BN  powder  and  those  calculated  in 
this  work 
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DISCUSSION 

The  experimental  results  and  theoretical  model  presented  above  indicate  that  we  are  capable  of 
experimentally  determining  and  theoretically  imderstanding  the  potential  of  various  powders  as  infrared 
obscurants.  Several  general,  but  important  observations  must  be  considered  about  the  work  conducted  so 
far.  The  data  presented  in  Fig.  4  verify  that  laboratory  measurements  on  a  small  scale  can  be  used  to 
predict  the  effectiveness  of  the  obscurants  when  deployed  in  a  large-scale  test.  The  ability  to  find  powders 
with  structured  resonances  throughout  the  IR  provides  the  confidence  that  we  can  choose  chemical 
characteristics  to  tailor  the  spectral  response  of  the  material  in  a  wide  variety  of  ways.  Finally,  the  succes 
of  the  theoretical  model  in  predicting  structures  and  cross  sections  will  provide  a  valuable  tool  in  screening 
future  candidates  and  their  appropriate  size  distributions.  In  this  section,  we  briefly  discuss  the  impact  of 
the  results,  some  recommendations  for  candidates,  and  directions  for  future  work  in  this  area. 

The  most  surprising  result  from  this  work  is  the  determination  of  the  cross  section  for  scattering  by 
brass.  Brass  was  chosen  as  a  baseline  system,  one  that  we  believed  was  not  extremely  useful  from  a  Navy 
point  of  view,  but  one  that  we  need  to  understand  because  it  is  widely  used  and  documented.  The  cross 
sections  we  measure  are  about  a  factor  of  16  lower  than  those  reported  (Edwards  1992).  However,  the 
ability  of  the  theoretical  model  to  predict  the  experimental  observations  within  10%  provides  convincing 
evidence  that  the  reported  values  are  incorrect.  The  reasons  for  this  error  are  not  yet  understood,  but  may 
relate  to  the  fact  that  the  original  data  were  taken  in  a  very  dense,  veiy  long  path  length  cloud,  which 
absorbed  >98%  of  the  incident  radiation.  Under  these  conditions,  several  possible  scenarios  emerge.  First, 
the  amount  of  radiation  on  the  IR  detector  might  be  close  to  the  noise  floor  of  the  system.  This  cannot  be 
proven  without  access  to  the  original  data.  Second,  multiple  scattering  effects  may  have  played  some  part 
in  the  experiment.  This  does  not  seem  to  be  likely  since  the  model  developed  in  our  work  indicates  that  this 
is  not  a  problem  in  our  ejq)eriment  where  the  mass  loadings  are  higher  (but  with  much  shorter  path 
lengths).  Third,  the  sample  we  obtained  may  not  be  identical  to  that  for  which  data  are  reported  in  the 
literature.  Because  the  obscuration  in  brass  is  based  solely  on  scattering,  the  particle  size  and  geometiy 
might  be  expected  to  play  a  large  role  in  determining  the  cross  section.  We  have  theoretically  looked  into 
some  of  these  problems,  and  differences  are  evident  for  spherical  particles  vs  ellipsoidal  particles  or  flat 
flakes.  However,  these  differences  are  not  large  enough  to  account  for  the  discrepancy.  Regardless  of  the 
reasons  for  the  differences,  it  is  clear  that  the  brass  flakes  are  not  as  effective  an  obscurant  as  is  popularl 
believed.  The  primary  reason  brass  was  chosen  over  graphite  by  the  Army  was  its  packing  density,  which 
is  about  three  times  that  of  graphite.  Since  the  cross  section  for  graphite  is  reported  to  be  between  1.4  and 
2.6  m^/gm,  and  we  measure  0.08  m^/gm  for  brass,  graphite  is  clearly  much  more  attractive.  We  plan  to 
acquire  various  graphite  powders  and  flakes  and  perform  an  identical  analysis  on  them  to  confirm  the 
reported  cross  section  measurements.  Should  these  experiments  agree  with  the  reported  values,  we  would 
recommend  the  graphite  as  an  alternative  to  brass  flakes.  However,  as  with  the  brass,  the  graphite  powders 
seem  to  rely  solely  on  scattering  for  obscuration,  with  no  obvious  resonances  in  the  3-5  |am  spectral  region. 

We  performed  similar  experimental  and  theoretical  studies  on  boron  nitride  powder,  which  was  one 
of  our  initial  candidates.  We  found  that  the  spectra  in  the  literature  (Doll  1994)  were  incorrect;  however, 
we  were  able  to  derive  good  optical  constants  to  model  its  absorption  and  scattering  characteristics.  Figure 
15  shows  the  agreement  between  theoiy  and  experiment.  While  boron  nitride  itself  does  not  have  strong 
resonances  in  the  3-5  pm  region,  it  does  exhibit  a  strong  resonance  at  7.4  pm  {oext  =  0.27  m^/gm)  which 
tails  into  the  10pm  region.  This  may,  therefore,  be  useful  in  the  long  wave  IR.  The  available  optical 
constants  cover  only  wavelengths  longer  than  5  pm,  but  the  Mie  calculation  includes  the  particle  size 
distribution  for  this  sample  provided  by  the  manufacturer.  Considering  the  complexity  of  the  problem,  the 
agreement  is  excellent.  This  also  confirms  that  our  model  is  capable  of  providing  useful  guidelines  for 
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predicting  accurate  cross  sections  for  both  scattering  and  absorption  features  through  the  entire  IR  region, 
provided  the  complex  index  of  rejB^ction  and  particle  size  distribution  are  available.  We  are  confident  that 
we  can  now  use  the  theory  to  prefilter  candidate  materials  and  provide  specific  guidelines  concerning  the 
optimum  particle  sizes  for  these  materials. 

The  data  shown  in  Fig.  5  for  boric  acid  reveal  several  absorption  peaks  in  the  regions  of  interest.  It 
should  be  noted  that  this  sample  was  not  analyzed  to  determine  its  size  distribution,  and  so  could  not 
accurately  be  modeled.  Although  the  spectra  appear  promising,  the  cross  sections  are  too  low  (typically  - 
0.008  m^/gm)  to  make  it  an  effective  obscurant.  However,  it  could  be  used  for  demonstration  purposes  or 
proof-of-principle  tests.  We  are  optimistic  that  this  is  not  the  best  system  available  in  this  region,  and  are 
looking  at  several  alternatives.  This  molecule  is  also  a  good  test  case  to  determine  our  ability  to  enhance 
the  cross  sections  by  increasing  the  scattering  part  of  the  extinction.  This  system  may  also  be  useful  as  a 
model  for  the  types  of  bonds  that  will  exhibit  strong  absorption  in  the  3-5  pm  and  8-12  pm  regions. 

Currently  we  are  screening  additional  materials  for  possible  application  to  this  problem.  Materials 
that  appear  promising  are  inorganic  amine  powders,  plastic  spheres,  graphitic  species,  and  some  minerals. 
Liquid  aerosols  may  also  be  promising  sources  as  possible  obscurants.  In  many  cases,  spectra  are  not 
available  for  these  species,  and  we  will  need  to  determine  these,  as  well.  It  does  appear  that  the 
methodology  and  approach  we  have  described  is  both  sound  and  practical  and  can  be  applied  to  a  nmnber 
of  problems.  We  also  propose  to  extend  this  technique  to  other  regions  of  the  spectrum,  including  the  near 
IR,  visible,  and  UV.  These  extensions  appear  to  be  straightforward  applications  of  the  current 
experimental  setup,  with  some  additional  optics  to  cover  the  expanded  portions  of  the  spectrum. 
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